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^ Abstract 

Recent theoretical work on the role of microscopic chaos in the dynamics and relaxation of 
^ many-body quantum systems has made several experimentally confirmed predictions about the 

+1j systems of interacting nuclear spins in solids, focusing, in particular, on the shapes of spin echo 

responses measured by nuclear magnetic resonance (NMR). These predictions were based on the 
idea that the transverse nuclear spin decays evolve in a manner governed at long times by the 

g 

1 slowest decaying eigenmode of the quantum system, analogous to a chaotic resonance in a classical 

^ system. The present paper extends the above investigations both theoretically and experimentally. 

I ^1 On the theoretical side, the notion of chaotic eigenmodes is used to make predictions about the 

CN relationships between the long-time oscillation phase of the nuclear free induction decay (FID) 

> 

and the amplitudes and phases of spin echoes. On the experimental side, the above predictions 

are tested for the nuclear spin decays of ^^F in CaF2 crystals and ^^^Xe in frozen xenon. Good 

CN agreement between the theory and the experiment is found. 
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I. INTRODUCTION 



The role of microscopic chaos in the observable behavior of macroscopic objects is a 
notoriously difficult elusive issud^. Exploring it requires joint analytical, numerical and 
experimental efforts. The focus of this paper is on a possible implication of microscopic chaos 
for free induction decays (FIDs) and spin echoes measured by nuclear magnetic resonance 
(NMR) in solids. 

Previously, one of us (B.V.F.) has arguecP^ that, as a consequence of microscopic chaos 
induced by generic non-linear interaction between nuclear spins, the long-time behavior of 
nuclear FIDs in solids has the universal long-time form 

F{t) = Ae-^^cos{iut - ifa), (1) 

where A, 7, u and cpa are constants. Depending on the microscopic Hamiltonian of interacting 
nuclear spins, the frequency u may be equal to zero. However, in the most common case of 
the magnetic dipole interaction, 00 has a finite value (see the discussion in Ref.tSf). Indeed, 
such behavior was observed as generic in the experimentaP^ and numericaP^ studies of 
quantum and classical spin systems. In a typical case, the constants 7 and u fall on the 
fastest natural microscopic time scale of the nuclear spin system, thereby precluding any 
explanation of the above behavior in terms of a damped harmonic oscillator — such an 
explanation would require the separation of time scales between the slow observable F{t) 
and much faster microscopic motion. 

The theoretical analysis of Ref.'^ predicted only the functional form ([l]) of the long-time 
FID behavior without predicting the parameters A, 7, u and A later paper^^ went 
further and predicted that different spin echoes initiated in the same system by perturbing 
the FID with almost any sequence of radio-frequency (rf ) pulses would have different initial 
behavior but then evolve to exhibit the long-time behavior ([T]) characterized by the same 
time constants 7 and 00. This prediction was confirmed experimentally in Refs.li^liil for 
hyperpolarized solid xenon and CaF2. 

The present paper explores the chaos-related notion that the long-time behavior of many- 
spin density matrices created in the course of the FID has a self-similar form accompanying 
exponentially decaying oscillations. In Section [IT], we show theoretically that for the spin 
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echoes initiated by perturbing the FID in the above long-time regime, all possible shapes 
of the echo responses are superpositions of two basics shapes with the relative weight of 
each shape determined by the phase of the FID oscillations at the time of the echo pulse. 
This two-shape decomposition is predicted to lead to a definite relationship between the 
long-time oscillation phase of the original FID and the long-time oscillation phases of the 



echo responses. In Section III, we verify the above predictions experimentally for CaF2, and 
include measurements on solid xenon in the appendix. 



II. THEORY 

The prediction of universal behavior of quantum spin systems in Refs.'^'^ was based on 
the conjecture that the long-time behavior ([T| is a manifestation of the slowest decaying 
chaotic eigenmode of the time-evolution operator, similar to a PoUicott-Ruelle resonanc d^ * ^^ *. 
Such eigenmodes control not just one observable quantity F{t) but also the evolutions of 
many-spin density matrices within the system 

Pki{t) = Po,fcze-(^+-)* + pi,,,e-(^— )* (2) 

where Pki{t) is the density matrix for any finite subsystem of the entire spin system, i.e. 
Pki{t) can be a one-spin density matrix, two-spin density matrix, or, in general, an n-spin 
density matrix, provided n is much smaller than the total number of spins in the system. As 
is often done in the NMR literature, Eq. ^ represents the leading correction to the infinite 
temperature density matrix pki = 1. The term pki = 1 does not contribute to the measured 
spin polarization. This high-temperature approximation should remain valid as long as 
the initial nuclear polarization is not too large, in the sense that the initial energy of the 
nuclear spin system with respect to the effective Hamiltonian of nuclear spin-spin interaction 
in the Larmor rotating reference fram^^ is close to the energy of the infinite temperature 
state. If the initial polarization is too large, then the system is expected to relax to a 
finite temperature equilibrium determined by its initial energy, in which case Eq.Q would 
represent the correction to the equilibrium density matrix for the final temperature. 

Equation ^ is the only connection between the theoretical analysis in this paper and 
the notion of microscopic chaos. Namely, the assumption of microscopic chaos justifies the 
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proposition that well-defined relaxational eigenmodes of the time-evolution operator of the 
entire system exist. In turn, the notion of an eigenmode of the time-evolution operator 
implies that the values of 7 and u do not depend on the order of the density matrix. What- 
ever the initial form of the n-spin density matrix, the long-time behavior would then be 
dominated by the slowest chaotic eigenmode of form ([2]) (among those compatible with the 
symmetry of the initial density matrix). The time-independent non-Hermitian form of po.fci 
for a given order of the density matrix, as well as the values of 7 and w, are determined by the 
microscopic Hamiltonian of the system. While the above connection to microscopic chaos is 
very indirect, we are not aware of any other framework justifying Eq.([2]). Our assumption of 
microscopic chaos is, in turn, motivated by the non-integrable character of nuclear dynamics 
governed by the nuclear spin-spin interaction Hamiltonian in the Larmor rotating reference 
fram^. 

The experimental evidence available so far is obtained from the total polarization of 
nuclear spins and, as such, indicates that the decay ([2]) is certainly present in the behavior of 
the one-spin density matrix, but not necessarily two-spin, three-spin and progressively higher- 
order spin density matrices. Higher-order density matrices are responsible for higher-order 
nuclear correlations (spin coherences in NMR language). The FID starts from a factorizable 
density matrix for the entire systeni^, meaning that the initial nuclear spin configuration is 
uncorrelated. Therefore, the expectation behind Eq. ([2]) is that the higher-order correlations 
first develop dynamically and then start decayin^^^, eventually approaching form (|2]). The 
effect of the echo pulse does not reverse but rather modifies the higher-order correlations. The 
predictions made below about the relationship between the shape of the echo response and the 
phase of the long-time FID beats at the time of the echo pulse are expected to be incorrect 
if the many-spin density matrices preceding the echo pulse do not exhibit the long-time 
behavior of form (|2| with the same parameters 7 and u independent of the order of the density 
matrix. On the other hand, the experimental confirmation of this relationship significantly 
strengthens the picture based on the notion of chaotic relaxation modes (PoUicott-Ruelle 
resonances) . 

We use the theoretical framework of Ref.^. The quantity called the "signal" is the total 
polarization of nuclear spins transverse to the external magnetic field. We consider the NMR 
response to the sequence of two rf pulses 
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90;-r-X (3) 

where the 90° pulse initiates the free induction decay (FID) and, ai^er the delay time r, 
pulse X "scrambles" the time evolution of the spin system. The FID between the two pulses 
is to be denoted by function F{t), and the signal at time (t — r) after the second pulse 
is to be characterized by the echo response function F(r, t). Time t is understood to be 
measured from the time of the first pulse. Most of the experimental tests of the FID-spin- 
echo relationships reported in Refs.'^^'^ used the solid echo pulse sequence characterized by 
X = 90«]. 

In the rest of this paper, we focus on the echo response F(r, t), which is initiated at a 
time r sufficiently long such that the FID function F{t) has already reached the asymptotic 
form ([T]). This regime was considered in Ref.l^, but there the main focus was on obtaining 
the envelope of the Hahn spin echo sequence 90°— r — 180°. — r for heteronuclear spin systems 
in an inhomogeneous magnetic field, in which case the echo can be monitored only at time 
T following the second pulse. In the present paper, we assume that the magnetic field is 
homogeneous, and thus that the echo response to the pulse sequence ^ can be monitored 
at any moment of time following the pulse X. We further assume that the shapes of the 
FIDs and echoes are determined by the dynamics of an isolated system of interacting nuclear 
spins in the Larmor rotating reference frame. The interaction Hamiltonian is assumed to be 
non-integrable, such as the case of the standard Hamiltonian of truncated magnetic-dipolar 
interaction!^^. 

When cu ^ 0, the long-time behavior of the density matrix ^ consists of the sum of 
the two Hermitian-conjugate terms po,kiG~^'^^^'^^'^ and p\ ki^~^"'~^'^^'^ ■ Each of these terms 
evolves in time in a self-similar way, in the sense that the evolution is controlled by the 
time-independent matrix po,fc« oi' p\ki^ while the time evolution of the entire density matrix 
Pki{t) is reduced to rescaling each of the above terms and changing their relative phase. As 
a result, we can also express the long-time behavior of the FID signal as the sum of two 
corresponding contributions 

F{r) = f{r) + nr), (4) 
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where 

fir) = ^ a e-(^+-)^ (5) 
and a is a complex-valued constant. Following pulse X, the new density matrix becomes 

Pki{r+) = e-(^+*")^f/xPo,M + e-^-^-^^^^Uxpiki^ (6) 
where Ux is the quantum operator representing the effect of pulse X. As a result, we obtain 

F{T,t) = f[r)f{t-r) + r{r)r(t-r) 



cos(a;r - v9„)Re/(t - r) + sin(u;r - V9a)lm/(t - r) , (7) 



where f{t — t) is the self-similar shape of the echo response associated with the first term 
in Eq. (|6]), and (fa is the complex phase of a. 

Equation ([T]) implies that one can experimentally measure any two echo responses F{ti, t) 
and F{T2,t), such that uj{t2 — ti) is not equal to a multiple of vr, then extract from these two 
responses functions Re/(t — r) and lTsif{t — r) and then, finally, predict F{t, t) for all other 
r. In fact, function Im/(t — r) can be measured directly by applying pulse X at a node of 
the FID, where cos(a;r — ipa) = 0, while Re/(t — r) can be measured by applying pulse X in 
the middle between two nodes, where sin(a;r — (fa) = 0. 

We can now elaborate on the long-time behavior of f{t — r) in order to relate the phases 
of the FID beats with the phases of the echo responses. The long-time behavior of f{t — r) 
is expected to be of the following fornP^ 

f(t -t) = 6ie-(^+'")(*-") + 62e-(^-'")(*~"), (8) 

where bi and 62 are two complex-valued constants which are not necessarily complex con- 
jugates of each other. (The only requiment here is that F(r, t) given by Eq. ([T]) is real.) 
Substitution of Eq. g into Eq. Q gives 

F{T,t) = ^ \a\ e-^* |C(r)| cos^ + y^cM], (9) 
where |C(r)| and v^cl^) are the amplitude and the complex phase of the function 

C{t) = ble-'"^^ + 626*^^""^'^"^ (10) 
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Both |C(r)| and <fc{T) should be independently accessible experimentally. The convenient 
representation of Eq.([Io| to test is 



|C(r)|' = + I62I' + 2|6i||62|cos(2a;r - 2^^ - Vb, - ^b,), (11) 

\C{t) \ cosipcir) = \bi\cos{ipa + v?6i) + |62|cos(2a;r - (pa - (Pbi), (12) 

where ipb^ and ipb^ are the complex phases of constants bi and 62, respectively. 

The FID function F{t) and the family of echoes F{t, t) can be measured experimentally 
choosing r and t — t large enough that the long-time regime ([T]) is reached for both F{t) 



and F{T,t). The test of Eqs. (11) and (12) can then be carried out in the following way: 

1) The parameters |a|, (pa, 7 and u are obtained from the FID asymptotics. 

2) The values of |C(r)| and (pciT) are obtained for each r by fitting the tails of the echo 
responses F(r, t) to Eq. (|9| as a function of t. 

3) Equations ( TTp!2 ) predict that both |C(r)p and |C(r)| cosv9c(r) consist of two terms: 



a r-independent constant and a r-dependent term oscillating with frequency 2uj. The ob- 
servation of this behavior as a function of r already constitutes a non-trivial test of the 
theory. 

4) The experimental curves for |C(r)p and |C(r)| cosipcir) are parameterized as follows: 

|C(r)p = fii + Di cos(2wr + 0i), (13) 

\C{t)\ cosifcir) = B2 + D2 cos(2u;r + ^2), (14) 

where the six parameters Bi, B2, Di, D2, (pi and (^2 should be directly accessible. The 
choice of the phases 0i and 02 are made such that the values of Di and D2 are positive. The 
theoretical formulas ( llp2 ) depend on four real- valued parameters: I62I5 fbi and iph2- 



Therefore, two further non-trivial tests are possible. 



Test 1: One obtains I62I = -D2, and then = y Bi — |62P- Prediction: Di = 2|6i||62|- 
Test 2: One obtains ipb2 = ~02 — fa, and then ipb-^ = —(pi — 2ipa — 'Pb2- Prediction: 

B2 = \bi\cOs{ipa + <fbi) = |6l|cos(02 " 0l)- 

As a final remark, we would like to mention that, in the case of monotonic long-time 
decays [w = in Eq.([T])], the treatment analogous to the one presented above predicts that 
the echo responses exhibit a self-similar shape with monotonic exponential long-time tails, 
which on a semi-logarithmic plots of the type of Fig. [T] would fall onto the same line. 



III. EXPERIMENT 



This section is focused on CaF2 where we were able to initiate echoes in the long-time 
regime of the FID. We also have the results from ^^^Xe in solid xenon, but in this material the 
latest echoes we were able to measure were not yet quite in the long-time regime of the FID. 
The are also additional theoretical complications related to the polycrystalline nature of the 
solid xenonl^that would make the predictions less rigorous, even if the echoes obtained were 
well into the long-time portion of the FID. As a consequence, the ^^^Xe data are included in 
an appendix. 

The FID and solid echoes of ^^F in CaF2 were acquired at room temperature in an 
external magnetic field of 2 T (^^F Larmor frequency 83.55 MHz). The CaF2 crystal used 
in our experiments was obtained from Optovac, Inc and is lightly doped with paramagnetic 
impurities (0.01% Ce) to reduce Ti to ~ 2 sec at 2 T. The CaF2 crystal was prepared with 
[100] axis along the long dimension of the cylinder. The crystal was then held with the 
magnetic field along the [001] direction. The data were acquired using the same Tecmag 
spectrometer using 2 /is square pulses with a receiver dead time of 13 /is. Using these 
parameters, 1000 transients were averaged with a repetition time of 10 seconds for each 
experiment to enhance the signal-to-noise ratio. 32 solid echoes were acquired, one every 
2.5 fis from 16 fis to 96 fis. At 60 fis, fits to the FID show that it has entered the long- 
time regime described by Eq. ([T]); therefore, echoes generated after time t = 60 /xs meet the 
criteria for testing the predictions made in Section [Tij The family of representative echoes is 
shown in Fig. [T] 

Determination of the amplitude |a|, decay coefficient 7, beat frequency u, and complex 
phase (fa of the FIDs in CaF2 was made by fits to Eq. ([T]) (see Table |l]). The amplitude |C(r)| 
and complex phase (fc of each echo signal were determined by fitting each solid echo signal 
to Eq. ([9]) with 7, to, \a\, and (fa fixed to the values obtained from the fit to the FID. 

In Fig. [2] we plot |C(r)p and |C(r)| cos^c" for the solid echoes of ^^F in CaF2. The solid 
lines are the fits to either Eq. (13) or Eq. (14), from which the parameters in Table [l] were 
obtained. Figure [2] shows the results for echoes initiated in both the early-time and the 
long-time regimes of the FID in order to illustrate the approach to the long-time behavior 
described by Eqs.( l3p4 ). 
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FIG. 1: (Color online). ^^F solid echoes in CaF2. Half of the 32 acquired solid echoes are 
shown, split between (a) and (b) for visual clarity. The echoes shown are acquired every 
5 fis from 16 /is to 96 /is. The first point of each echo is indicated by a solid circle to guide 
the eye. 

Tests 1 and 2 formulated at the end of Section [IT] are then carried out. The predicted and 
the measured values of parameters Di and B2 are listed in Table |TTj We find that in each 
test the predicted and the measured values agree with each other within the experimental 
uncertainties. 

Finally, Section [IT] contains a more general prediction: that all possible shapes of echo 
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FID Parameters 


Solid Echo Parameters 


\a\ = 9, 380, 000 ±45, 000 
<fa = 1.921 ±0.006 
7 = 0.0414 ± .0008 /xs-i 
a; = 0.120 ±.007 ^s-^ 


Si = 18. 3 ±0.2 ^2 = 0.02 ±0.3 
Di = 2.1 ±0.3 D2 = 4.2 ±0.4 
01 = -1.5 ±0.1 02 = 0.4 ±0.1 



TABLE I: Long-time FID and echo fit parameters for CaF2. The FID parameters are 
obtained by fitting the FID to Eq. ([!]). The echo parameters are obtained by fitting the 



amplitudes and phases of the measured solid echoes to Eqs.(13) and (14). 



Test 1 


Test 2 


Mispredicted — 2.3 ± 0.3 
-Ml,measured 2.1 ± 0.3 


B2,predicted = "0.01 ± 0.07 
B2,measured = 0-00 ± 0.3 



TABLE II: Two tests of the theoretical predictions formulated at the end of Section 



II 



responses including both the initial and the long-time behavior can be obtained from two 
basic functions Re/(t — r) and Im/(t — r) as described by Eq. ([T]). To determine these 
functions, we chose two measured echo responses F{Ti,t) and F(r2,t) initiated at times 
Ti = 81.0 /is and T2 = 93.5 fxs, respectively. According to Eq. ([T]) 



F(ri, t' + n) = AiRe/(t') ± 5ilm/(t') 
F{t2, t' ± T2) = A2Ref{t') ± B2lmJ{t') 



(15) 
(16) 



where in Eq.(15) t' = t — ri, Ai = F{ti) and Bi = F(ri) tan(a;ri — </?„), while in Eq.(15) 
t' = t — T2, Ai = F{t2) and Bi = F{t2) taji{ujT2 — ^a)- Here Fiji) and F{t2) are the measured 
values of the FID at times ti and r2, respectively. Now we express functions Re/(t') and 
Im/(t') in terms of the measured functions F(ri,t' ± Ti) and F{T2,t' + T2) by solving the 
system of linear equations ( IsflG ) and then substitute Re/(t — r) and Im/(t — r) back to 
Eq.Q to predict other echo responses initiated in the long-time FID regime. The functions 
Refit') and Im/(t') are plotted in Fig. [sj Since for solid echoes F{t,t) = F{t), Eq.^ 
implies that Re/(0) = 1 and Im/(0) = 0; however, due to the finite (13 /is) recovery time of 
our spectrometer, we were unable to measure these functions back to time t' = 0. 
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FIG. 2: (Color online) (a) Amplitudes and (b) phases of ^^F solid echoes in CaF2 as a 
function of inter-pulse delay time r. Open circles represent echoes generated in the 
early-time of the FID, while solid circles represent echoes generated in the long-time of the 



FID. The solid line (red) is the best fit of the long-time data to Eq. (13) for (a) or to 



Eq. (14) for (b). 



In Fig. |4| we show several of the measured echoes in CaF2. Each echo F{t, t' + r) in Fig. |4| 
has been multiplied by a factor e'^'^ to correct for the exponential decay of the FID at their 
respective initial time values. Labeled in the plot legend are the two echoes F(ri,t' + ri) 
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FIG. 3: (Color online) The shape functions Re/(t') and Im/(t') obtained from the linear 



system of equations (15 16). 



and F{T2,t' + T2) used in determining the shape functions Ref(t') and lmf{t'). In Fig. [sj the 
measured echo responses initiated at r = 61,66,76,83.5,86, and 91 fis are compared with 
the predicted ones obtained by substituting Re/(t') and Im/(t') into Eq.([7|. We observe in 
Fig. [5] that the agreement of the long-time behavior between the predicted and the measured 
echo responses is very good for all echoes, which is consistent with the results presented in 
Table [n] and in Fig. |2} The initial behavior of the early echo responses (r = 61 /is, 66 /is, 
and 76 /is) exhibits some discrepancies between the predicted and the measured behavior. 
However, it is clear that the predicted behavior still captures the evolution of the measured 
echo shapes in a satisfactory way. In the later echoes (r = 83.51 /is, 86 /is, and 91 fis), this 
initial discrepancy no longer appears and the entire echo shape is found to agree with the 
predicted shape. 

A possible reason for the above initial discrepancy in the early echoes may be the presence 
of the chaotic eigenmodes that decay faster than the slowest mode controlling the long-time 
FID behavior ([T]), but still not fast enough to completely disappear by time r when the echo 
is initiated. The existence of a well-separated second slowest eigenmode was demonstrated 
by the recent experiment of Meier et al.l^ on a CaF2 crystal for the same orientation with 
respect to the magnetic field. This second eigenmode disappears below the noise level on the 
timescale of 60 /xs. These additional eigenmodes are probably more pronounced in the higher- 
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t' (us) 

FIG. 4: (Color online) ^^F solid echoes in CaF2 labeled by their interpulse delay times r. 
The quantity plotted is e^'^F{T,t' + r). The 81.0 /is and 93.5 /is echoes represent the echoes 
used to obtain the shape functions Ref(t') and lmf(t'). 



order spin correlations, because these correlations develop in the course of the FID evolution 
only after an initial time delay with respect to the beginning of the FID^^. Therefore, the 
behavior of many-spin density matrices should approach the long-time form ^ also with 
some delay with respect to the time when the FID starts exhibiting the universal long-time 
form ([T]). The same additional eigenmode may be controlling the approach of |C(r)p to the 



predicted asymptotic behavior in Fig. 2a 



IV. SUMMARY AND CONCLUSIONS 



In this paper, we investigated the properties of spin echoes initiated in the regime where 
the nuclear FID has reached the universal exponentially damped oscillatory behavior. Using 
the theoretical framework motivated by the notion of microscopic chaos, we predicted how 
the shapes of the echo responses depend on the phase of the FID oscillations at the time of 
the echo pulse, and, in particular, obtained the phase relationships between the long-time 
oscillation of the FIDs and the echoes. We have further conducted several experimental tests 
of the above predictions for FIDs and solid echoes in CaF2 and solid xenon, and obtained 
results in good overall agreement with the theoretical expectations. The long-time phase 
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FIG. 5: (Color online) ^^F solid echoes in CaF2 (solid lines). The red (broken) lines show 
the theoretical echo shapes obtained using the Refit') and lmf(t') shape functions. The 
quantity plotted is e^'^Fij^t' + r). The values of r are indicated in the plot legend. 



relationships between the FID and the echoes were confirmed particularly well. This good 
agreement amounts to an indication that the long-time behavior of the higher-order spin 
density matrices has the form given by Eq. (g with the same values of 7 and uj as the 
original FID. Such a behavior is expected for a relaxational eigenmode of the time-evolution 
operator in a chaotic system. 

While the fundamental difficulties in defining the notion of microscopic chaos still remain, 
the present paper demonstrates that the approach of Refs.'^'^ based on making parallels with 
relaxational eigenmodes in classical chaotic systems continues to generate successful quan- 
titative predictions. These predictions were made in a regime not accessible by controllable 
first principles calculations. We are not aware of any other approach that would reproduce 
the same predictions under conditions that the quantities of interest (nuclear spin decays) 
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evolve on the fastest microscopic time scale of the system. 
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Appendix A: Measurements in Solid Xenon 

We have carried out measurements on ^^^Xe in solid xenon similar to those reported 



in Sec. Ill on CaF2, but due to experimental constraints described below, we were not 
able to generate solid echoes in the long-time regime of the FID [i.e. where the FID is 
well-described by Eq. ([T])] that themselves had enough signal-to-noise ratio to be measured 
into their respective long-time regimes. However, as the echoes acquired with the longest 
interpulse delays show an approach to the predicted forms, we present the data even though, 
strictly speaking, they were not generated in a regime wherein the predictions hold. 

For this appendix we use the FID and solid echoes which are presented in Ref.'^. Poly- 
crystalline xenon samples were produced in a magnetic field of 2 T (^^^Xe Larmor frequency 
24.56 MHz) using the methods described in Ref.'^. Ten solid echoes were acquired approx- 
imately 0.2 ms apart from 0.4 ms to 2.5 ms. A fit of Eq. ([T| to the FID show that it does 
not enter the long-time regime described by Eq. ([T| until after t = 2.5 ms. As the latest 
echo was acquired with an interpulse delay time of 2.5 ms, no echoes were acquired in the 
long-time regime of the FID. The xenon FID and solid echoes are shown in Fig. [6] 

The amplitude \a\, decay coefficient 7, beat frequency u, and complex phase ipa of the 

129^ 



'Xe FID were determined by fits to Eq. ([T]) (see Table HI). The amplitude |C(r)| and 
complex phase (pc of each echo signal were determined by fitting each solid echo signal to 
Eq. ^ with 7, u, \a\, and (fa fixed to the values obtained from the fit to the FID. 

In Fig. ul we plot |C(r)p and | C(r )| cos 0c for the measured solid echoes. The solid 



lines are the fits to either Eq. (13) or Eq. (14), from which the parameters in Table HI 



were obtained. Figure [7] shows the results for echoes initiated in both the early-time and 
the long-time regimes of the FID. We observe that the latest echoes begin to approach the 



behavior predicted by Eqs. (13) and (14). These echoes are labeled "late-time echoes" even 
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FIG. 6: (Color online). ^^^Xe solid echoes in solid xenon. The latest 4 echoes acquired are 
shown. The first point of each echo is indicated by a solid circle to guide the eye. 



thought they are not actually in the late-time region as is the case in the CaF2. 

Tests 1 and 2 formulated at the end of Section [IT] are then carried out. The predicted and 



the measured values of parameters Di and B2 are listed in Table IV In each test, the pre- 
dicted and the measured values agree with each other within the experimental uncertainties. 

Finally, we compare the obtained echoes with the predicted shape functions. We first 
obtain the shape functions Re/(t') and lmf{t') as described in Sec. 



Ill 



In Fig. 



8a 



we show 



the measured echoes in ^^^Xe. Each echo F{t, t' + r) in the figure has been multiplied by 
a factor e''"^ to correct for the exponential decay of the FID at their respective initial time 



FID Parameters 


Solid Echo Parameters 


\a\ = 151, 800 ±1,600 
ifa = -1.254 ±0.006 
7 = 1.251 ±0.005 ms'i 
uj = 2.10 ±0.01 ms-^ 


Si = 260 ±30 B2 = 0.1±3 
Di = 290 ±40 = 12 ± 3 
01 = 0.71 ±0.05 (/.2 = 2.4 ±0.5 



TABLE III: Long-time FID and echo fit parameters for solid xenon. The FID parameters 
are obtained by fitting the FID to Eq. ([T]). The echo parameters are obtained by fitting the 



amplitudes and phases of the measured solid echoes to Eqs.(13) and (14). 
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o Early-time echoes 
• Late-time echoes 
Best ft to tlneory 




(b) 



FIG. 7: (Color online) (a) Amplitudes and (b) phases of ^^^Xe solid echoes in solid xenon 
as a function of inter-pulse delay time r. Open circles represent echoes generated in the 
early-time of the FID, while solid circles represent echoes generated in the long-time of the 



FID. The solid line (red) is the best fit of the long-time data to Eq. ( 13 ) for (a) or to 



Eq. (14) for (b). 



Test 1 


Test 2 


Dl,predicted = 250 it 60 
D 1 ^measured — 290 it 40 


B2, predicted = —1.0 it 0.3 
B2,measured — it 3 



TABLE IV: Two tests of the theoretical predictions formulated at the end of Section 



II 
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FIG. 8: (Color online) (a) ^^^Xe solid echoes in solid xenon labeled by their interpulse 
delay times r. The quantity plotted is e'^'^F{T,t' + r). The 1.7 ms and 2.5 ms echoes 
represent the echoes used to obtain the shape functions Ref{t') and Im/(t'). (b) The 



functions Re/(t') and Im/(t') obtained from the linear system of equations (15 



16). 



values. The two echoes F{Ti,t' + ri) and F{T2,t' + T2) used for determining the functions 
Re/(t') and Im/(t') are shown in Fig. 8a, The functions Re/(t') and Im/(t') are shown in 
Fig.|8b} 

In Fig. |9l the remaining two late-time echoes are compared with the predicted ones 



obtained by substituting Re/(t') and Im/(t') into Eq.([7p. We emphasize that the theoretical 
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prediction of the echo shapes is expected to hold only for echoes initiated after 2.5 ms in 
this material. Therefore, the discrepancy in the early parts of the echo shapes is expected 



for the reasons discussed in Sec. Ill in relation to the early-time echo shapes in CaF2. 




2 3 4 
Time (p-s) 

FIG. 9: (Color online) ^^^Xe solid echoes in solid xenon (solid lines). The red (broken) lines 
show the theoretical echo shapes obtained using the Re/(t') and Im/(t') shape functions. 
The quantity plotted is e^'^F{T,t' + r). The values of r are indicated in the plot legend. 



We finally mention a possible additional complication in solid xenon associated with 
the fact that our xenon samples are not single crystals but rather polycrystallites. Our 
very recent theoretical studj^ indicates that the observed long-time FID behavior of the 
polycrystalline solid xenon probably represents a typical long-time behavior of the individual 
crystallites contributing to the polycrystallite average, but the true asymptotic FID behavior 
of the entire polycrystallite is expected to appear only at times beyond the range accessible 
in our experiments. This asymptotic FID behavior should be controlled by the small fraction 
of the constituent crystallites with the slowest exponential decay constants 7. 
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